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Abstract 

The laser driven acceleration of ions is considered a promising candidate for an ion source for hadron therapy of oncological 
diseases. Though proton and carbon ion sources are conventionally used for therapy, other light ions can also be utilized. 
Whereas carbon ions require 400 MeV per nucleon to reach the same penetration depth as 250 MeV protons, helium ions 
require only 250 MeV per nucleon, which is the lowest energy per nucleon among the light ions. This fact along with the larger 
biological damage to cancer cells achieved by helium ions, than that by protons, makes this species an interesting candidate for 
the laser driven ion source. Two mechanisms (Magnetic Vortex Acceleration and hole-boring Radiation Pressure Acceleration) 
of PW-class laser driven ion acceleration from liquid and gaseous helium targets are studied with the goal of producing 250 
MeV per nucleon helium ion beams that meet the hadron therapy requirements. We show that He® ions, having almost the 
same penetration depth as He^ with the same energy per nucleon, require less laser power to be accelerated to the required 
energy for the hadron therapy. 

PACS numbers: 52.25.0s, 52.38.Kd, 52.27.Ny 
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I. INTRODUCTION. 


Beams of high energy particles, initially produced for the study of fundamental processes of particle physics, found 
applications in different areas of science and technology. In particular, the medical applications of such beams gave 
rise to the field of radiation therapy [ij. Hadron therapy, the part of radiation therapy that uses ion beams has a 
number of advantages, such as a dose profile with a high dose deposition near the beam stopping point (Bragg peak). 
It provides a conformal dose deposition treatment, in which the radiation dose is matched to the tumor, while the 
radiation induced damage to the normal tissue is limited. It is known that the light ions (such as protons) allow 
for a more precise dose delivery in the longitudinal direction, but experience more side scattering than heavier ions 
(such as carbon). Heavier ions, additionally, generate higher radiation induced damage (RBE - Relative Biological 
Effectiveness). This damage is hard to compensate by the cell repair mechanisms [^, which is critically important for 
the treatment of radioresistant tumors. Heavier ions, however, generate a tail in the energy transfer distribution near 
the beam stopping point mainly due to the fragmentation of ions in the course of interaction with the atom nuclei in 
the tumor. The fragmentation may lead to the unwanted irradiation of surrounding healthy tissue. This effect can be 
mitigated by using the Intensity-Modulated Radiation Therapy. Mostly protons and carbon ions are used at hadron 
therapy facilities, however different species of ions have also been tested. One can consider optimizing ion species to 
maximize the RBE for the necessary accuracy of the dose delivery. 

Moreover, hadron therapy may be optimized not only by using different ion species but also by using different 
ion sources. Presently all hadron therapy facilities in operation are based on the conventional accelerator technology 
making use of accelerators and complex ion beam transport systems (gantry) for beam delivery in the treatment 
room [I|. It is widely discussed whether a different source of ion beams can be utilized to reduce the cost and size 
of these facilities. One of the possibilities is the use of high power laser systems Q, which, in principle, exhibit 
lower construction and operation costs, and moreover would require smaller beam transport system, leading to a 
more compact facilities. Due to the rapid evolution of laser technology it became possible to use such systems for 
efficient acceleration of charged particles The theoretical and computer modeling results show that PW-class 

lasers, which are already in operation @ and are being constructed are able to produce ions with the energy high 
enough to meet the requirements for hadron therapy [7|. Experimentally proton beams with the maximum energy of 
45 MeV for short (tens of femtoseconds) [^, Q and 160 MeV for long (hundreds of femtoseconds) laser pulses were 
observed. The proof-of-principle experiments testing the ability of the laser produced proton beams to damage cancer 
cells were recently carried out m- However the generation of a laser based ion source that would meet the maximum 
energy, energy spread, number of particles per bunch, and repetition rate requirements [mil simultaneously is still 
a challenge both theoretically and experimentally. 

In this paper we study different mechanisms of generating a laser driven source of helium ions. Usually either protons 
or carbon ions are considered for medical applications of laser acceleration [mil. Here we consider an intermediate 
ion species that may be more beneficial for therapy due to higher biological damage than that of protons and lower 
laser energy input than that for carbon ions. There are some indications that helium ions demonstrate higher biological 
damage and good dose localization with the lower dose deposition at the beam entry 0. The penetration depth of 
an ion beam is determined by its initial energy. For a d-shaped initial ion e nerg y distribution, /(0,£’) = d{£ — £a), 
the location of the Linear Energy Transfer (LET) maximum is Xmax = £aM [13j where rj = A'Ke^niaZuZ'^lme: here 
Zq, rua are the charge state and mass of the ion; is the charge state of the target; mg and e are electron mass and 
charge; £a is the initial ion energy. This relation allows us to estimate how much energy per nucleon is needed for an 
ion beam to have the same penetration depth as a proton beam with initial energy of £p: 

{£M = £p , ( 1 ) 

1 /2 

here Aa is the ion atomic number. For all stable isotopes of elements from helium to carbon the factor Za/AJ ranges 
from 1 to 1.73, being the smallest for He^ and the largest for carbon. The isotope He^ has Z}jg3 j= 1.15. This 
makes He^ and He^ ions interesting candidates for hadron therapy, since they have the lowest energy per nucleon 
requirement among the ions heavier than protons, and deliver higher biological damage. 

In order to be practical the laser driven ion source should be able to deliver helium ions with the energy up to 250 
MeV/u with a short pulse PW-class laser system and the target should be replenishable to allow for the high repetition 
rate. In the case of Helium it is reasonable to consider two types of targets: (i) liquid and (ii) gaseous. The first target 
can be realized as a liquid jet or a wall of a liquid helium bubble. The gas target may be realized as either a very 
narrow high density gas jet or a shock launched inside a low density gas jet. In principle, the available laser power 
and the properties of the target make it possible to choose the regime of interaction that would allow to optimize 
the parameters of the produced ion beam. There are several basic laser ion acceleration mechanisms: (i) Target 
Normal Sheath Acceleration (TNSA) [2, (ii) Coulomb Explosion (CE) 0 , (iii) Radiation Pressure Acceleration 
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(RPA) [13, (iv) Magnetic Vortex Acceleration (MVA) MM, and (v) the Shock Wave Acceleration (SWA) [^ . 
There are also several mechanisms that throngh either modification or combination of some of the basic regimes 
enhance the maximnm ion energy, nnmber of accelerated ions, or improve their spectrum. For example, the Burn-out- 
Afterburner (BOA) [13, which employes an enhanced TNSA, and Directed Coulomb Explosion (DCE) [13, which is 
the combination of RPA and CE, are such composite mechanisms. Also the use of composite targets (consisting of 
low density and high density parts) was proposed in a number of papers to either inject the ions into accelerating 
fields or to enhance the interaction of the laser pulse with the high density part of the target [13-113 ■ Most of the 
above mentioned mechanisms were shown theoretically [3 to produce high energy ion beams from either ultra-thin 
(tens or hundreds of nanometers) solid density targets or employing long pulse lasers, which can not be achieved in 
the case of liquid and gaseous helium targets and short laser pulses. Thus we consider two mechanisms of laser ion 
acceleration relevant to these targets and able to provide necessary peak energy. In the case of a very thin liquid 
target it is Radiation Pressure Acceleration (RPA) [13 US (13, ^ S®'® target it is Magnetic Vortex 

Acceleration (MVA) [H, [l3,113 ■ 


II. LIQUID HELIUM TARGET. RADIATION PRESSUREACCELERATION. 

The interaction of a laser pulse with a liquid helium target can be viewed as "hole-boring" [13, which is RPA of 
a density modulation inside an extended target. In principle, the energy gain in this regime can be estimated by 
assuming that there exists a stationary solution {dfi/dt = 0, where /3 is the velocity of the laser plasma interface). 
The position of the laser plasma interface is then determined by the balance between the ion momentum flux and the 
flux of the EM wave momentum. This gives [13 for the velocity of this interface: 

+ -B%(1+ Pi) 

- w^> -■ 

where Be = (E^/27rne?7i//eC^) Pg is the group velocity of the laser pulse, and rie is the electron density. Here we 
took into account the fact that a tightly focused laser pulse has a group velocity smaller than the vacuum light speed, 
which ultimately limits the maximum achievable ion energy [13, [13- For ^ 1 PW laser pulse focused down to 1.5/im 
focal spot at helium plasma with the density of 28ncr, Eq. m, where the laser pulse EM intensity is averaged over 
time and space, gives an energy gain of Eue — 0.6 GeV, which is ^ 150 MeV per nucleon. Here ricr = w^mg/dTre^ is 
the critical plasma density and w is the laser frequency. The scaling for the energy gain with the increase of the laser 
power from Eq. © is Eue ^ 

In what follows we present the results of 2D PIC simulations of a PW-class laser interaction with a He target with 
thickness ^ 6A, where A is the laser wavelength, and electron density of 28ncr for He'’’ and ISucr for He^, which 
corresponds to the liquid Helium density near the boiling point (4.7 K). The simulations are performed using code 
REMP M- The simulations are set up as follows: the simulation box is 20A x lOA, dt=0.0025A, and dx=dy=0.005A, 
chosen in order to resolve the skin depth. The laser is introduced at the left boundary with Super Gaussian longitudinal 
and transverse profiles, exp j — jii/w)^ — [2{x — t)/T]®}, duration of r = 30 fs, and spot size of w = 2A. The utilization 
of such laser profile was shown [22] to be beneficial for enhancing the maximum ion energy. The laser pulse is focused 
8A from the left border with f/D=2, giving the laser spot size at focus of about 1.5A. The laser pulse bores a hole 
in the liquid target, accelerating helium ions at the laser-plasma interface, which is a characteristic feature of the 
hole-boring RPA mechanism. The spectra of He® and He^ are exponentially decaying with a cutoff energy of 300 
MeV for 1 PW laser pulse and ^ 700 MeV for 2 PW laser pulse (see Fig. [T]). Only in the 2 PW case the helium ion 
energy partly overlaps with the required for hadron therapy energy range of (240 MeV,860 MeV) for He® and (280 
MeV, 1 GeV) for He^. However the fast decay of the spectra leads to a very small number of particles near the cutoff 
energy: 3 x 10^ in the energy interval {O.OSfj^e, l-02£^//e}, with Effg = 600 MeV. The utilization of thin liquid targets 
will require a high level of control over the laser prepulse to prevent the significant modification of the target before 
the arrival of the main pulse. In some cases the existence of the prepulse can boost the acceleration mechanism by 
providing more favorable conditions for the laser pulse coupling to the target, however high level of control over the 
prepulse is mandatory in this case [13 ■ 


III. GASEOUS TARGET. MAGNETIC VORTEX ACCELERATION. 

In the MVA regime, during the laser pulse propagation in near critical density (NGD) plasma the ponderomotive 
force of the laser expels electrons and ions in the transverse direction. The time of ion response is much longer than 
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FIG. 1. The spectra of helium-4 (blue curves) and helium-3 (red curves) ions accelerated from a 6A liquid target by 1 PW 
(dashed curves) and 2 PW (solid curves) laser pulses in the forward direction, i.e., pHe > 0. 


that of electrons, which leads to a formation of a positively charged region just after the front of the laser pulse. Thus 
the electrons move under the action of the laser pulse EM field and the field of a remaining ion core. This results 
in the formation of a channel in electron density with high density walls. After the laser pulse passes, the ions also 
expand in the transverse direction forming a channel in ion density. The propagation of an intense laser pulse in the 
relativistically underdense plasma can be approximated by the propagation of an EM wave in a waveguide 0- The 
radius of the self-generated channel can be estimated from balancing the electron energy gain from the laser pulse EM 

1 7*2 

field and from field generated by the positive charge of ions in the channel: Rch = {achncr/ne) ' A/tt. The relation 
between the amplitude of laser pulse vector-potential in the channel, ach, and the laser pulse power, P, is 


\ K Pc ‘^cr J 


(3) 


where Pc = 2m\(^ je^ = 17 GW is a characteristic power for relativistic self focusing [s^. The factor K = 0.074 
comes from the integration over transverse coordinates and the duration of the laser pulse. The radius of the channel 
in terms of the laser pulse power is 


Rch 


^\ne) \KPc) 


(4) 


In the MVA regime the acceleration is achieved by an electric field at the back of the target, which is generated 
by, and is of order of, the magnetic field {Bch) of electrons, accelerated by the laser pulse in plasma in the forward 
direction, as they exit the target. This EM field associated with the electron current leaving the target expands along 
the back surface of the target, however the magnetic field flux is conserved. The maximum energy gain of a charged 
particle in such an expanding field is Eue ^ ^ZneBchRch- 

Let us assume that the electron current consists of all the electrons that initially were in the volume which is 

the volume of a cavity left in the wake of a laser pulse, ensuring that the fields of electrons and ions are compensated 
outside the cavity and the strong fields exist only inside it. As the electrons move through the charged ion background 
they experience a plasma lensing effect (^ . i.e., they are pinched towards the central axis and the radius of the electron 
beams is determined from the balance of the transverse electric field of the ion column, Ei = 2'KencRb, which in the 
reference frame of an e-beam is E[ = 'jeEi, and the self field of an electron beam, which is equal to = 2'Ken'^Rb. 
From the condition E[ = A' we obtain Rb = Rch he- Here 7 e is the Lorentz factor of the bulk of electrons accelerated 
forward. In the regime of a laser pulse interaction with a NCD plasma the electrons are continuously injected in 
the cavity behind the laser pulse leading to a formation of an electron current, which is dominated by low energy 
electrons. 

The characteristic energy of these electrons is the injection energy, which can be obtained from the condition on 
the electron velocity to be equal to the group velocity of an EM pulse propagating in a waveguide of radius Rch'- 
7 e = (v^/l.Sd) {2P/KPc)^^^ {ucrlne)^^^■ The magnetic field, generated by these electrons, is Bch{Rb) = 2'KencRchle 
aX r = Rb- Then the energy gain of an ion in the electric field at the back of the target is 

Ehc — TficC 27r Zpic {p-cl^cr^ {Rch/^^ j (5) 
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He'' 

H^ 

P=500 TW 

8 X 10' 

380 MeV 

0 

190 MeV 

3 X 10* 

570 MeV 

0 

230 MeV 

1 X 10* 

400 MeV 

3 X 10' 

290 MeV 

P=1 PW 

6 X 10“ 

520 MeV 

2 X 10“ 
380 MeV 

9 X 10* 
670 MeV 

4 X 10* 
450 MeV 

8 X 10* 
830 MeV 

3 X 10* 

530 MeV 

P=2 PW 

1.6 X lO** 
770 MeV 

1.5 X lO** 
550 MeV 

2.6 X lO'' 
930 MeV 

1.2 X 10'' 
650 MeV 

1.5 X 10“ 
1100 MeV 

8 X 10* 
880 MeV 


TABLE I. The number of He^ and He"^ ions in the energy interval {0.98£max, l-02£max}, where £max = 250AHe MeV, and the 
maximum helium ion energy per nucleon vs laser power and density of the target. 



FIG. 2. The distribution of He"^ ion density in the course of a 1 PW laser pulse interaction with a GOA NCD He^ target at 
t — 130 X 2-k jijj. 


which gives an upper estimate for the ion energy gain in the MVA regime. For example, for a 1 PW laser pulse and 
rie = Sricr ^He ~ 3 GeV. From Eq. ([S]) we see that the ion energy scales with the laser pulse power as . The 
scaling was derived assuming the optimal match between the laser pulse and the plasma, i.e., the target length is 
equal to the laser pulse depletion length. We should mention here that MVA regime is relevant for ion acceleration 
from targets with densities being not far from critical (T^ . For lower densities, laser pulse filamentation and Langmuir 
wave generation prevent efficient channel generation. At higher densities other mechanisms of laser ion acceleration 
come into play and the target thickness, which is equal to the depletion length, is of the order the pulse length or 
smaller, which is outside the regime of MVA applicability. 

In what follows we present the results of 2D PIC simulations (sij of a PW-class laser interaction with He^ and 
He^ targets with thickness (30 - 100)A and electron density of (l-3)ncr, which corresponds to the density of a high 
pressure Helium gas jet. The simulations are set up as follows: the simulation box is 150A x lOOA, dt=0.05A, and 
dx=dy=0.1A. The laser is introduced at the left boundary with Super Gaussian longitudinal and transverse profiles, 
exp{ —(y/w)® — [2(x — tj/r]®}, duration of 30 fs, and spot size of w = 4A. The laser pulse is focused 16A from the 
left border with f/D=2, giving the laser spot size at focus of about 1.5A. The distribution of helium ions during the 




FIG. 3. a) The spectrum of helium-3 (red) from 240 MeV to 860 MeV and helium-4 (blue) ions from 280 MeV to 1 GeV 
accelerated in the forward direction; b) The number of helium-3 (red) and helium-4 (blue) ions per energy interval, 0.98fue < 
fife < 1.02£He- The laser pulse power is 1 PW, the target thickness is 60 A. 
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interaction is shown in Fig. [5J A thin ion filament formed along the laser propagation axis is clearly seen, as the ions 
from it are accelerated by the charge separation electric field formed at the back of the target. The spectra of He^ 
and He^ ions are shown in Fig. [3^. for the interval of energies relevant to hadron therapy: 280 MeV< <1 GeV 
and 240 MeV< fweS < 860 MeV. Such energy selection can be achieved by a system based on permanent dipoles as 
shown in Ref. [33 ■ The slow decay of the spectrum in the interval relevant to therapy needs makes it possible to 
utilize different energy beamlets corresponding to different penetration depths simultaneously, thus greatly increasing 
the effectiveness of the laser driven ion source. The number of ions per energy interval {0.98£He, 1.02fi/e} is shown 
in Fig. IHJr. We see that it demonstrates a very slow decay with a difference between the maximum and minimum 
number less than two times. 

We summarize the results of PIC simulations in Table I. The maximum energy per nucleon varies from 400 MeV 
to 1.1 GeV per nucleon for He^ and from 190 MeV to 880 MeV per nucleon for He"*^. He^ ions have higher maximum 
energy per nucleon for all considered density and laser power values, yielding a higher number of particles in the energy 
interval, {0.98£max, ^■9‘^£max}- According to simulation results, 500 TW laser power can be considered a threshold of 
therapeutic beam production with a number of particles per bunch of about 10® in the case of He®. Further increase 
of the laser power up to 2 PW shows an increased production of He ions with a number of particles per bunch reaching 
2 X 10^. 


IV. CONCLUSIONS. 

We considered a laser driven helium ion source for the hadron therapy. Two types of targets were studied: a liquid 
target (accelerated using hole-boring RPA) and a gaseous target (accelerated using RPA). We found that the liquid 
helium target requires much higher laser power to produce ion beams that may be of interest for therapy than the 
gaseous target, namely, 2 PW vs 500 TW, and produces significantly less particles per energy bin, namely, 10^ vs 
10®. Moreover the spectrum of the helium ions produced through the MVA regime allows for the multiple energy bin 
extractions to simultaneously treat the parts of tumor corresponding to different penetration depths. The number of 
particles per energy bin is about 10®, which combined with the 1 Hz laser repetition rate meets the hadron therapy 
requirements regarding the total dose and the procedure duration. We also investigated the feasibility of using He® 
isotope for the laser driven ion source. He®, having almost the same penetration depth as He^ with the same energy 
per nucleon, requires less laser power to be accelerated to the required energy for hadron therapy. This He® advantage 
should be weighted against the relative biological damage of He® and He'’’ ions, but this optimization lies outside the 
scope of this paper. 
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